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1. Introduction 

During the last two decades Supersymmetric Quantum Mechanics became an important 
and popular tool in establishing connections between a variety of isospectral quantum 
systems |]1| . It gives new insight into the problem of spectral equivalence of Hamiltonians, 
which historically was developed as Factorization Method of Schrodinger in Quantum 
Mechanics and as Darboux-Crum transformations in Mathematical Physics |^. 

In standard 1-dim SUSY QM |Q two almost isospectral Hamiltonians Hi and H2 are 
combined into a single (matrix) Schrodinger operator in a Hilbert superspace. The two 
SUSY partners of a Superhamiltonian are intertwined by supercharges Q^ 

HiQ+ = Q+H2] Q-Hi = H2Q-; Q- = {Q+Y (1) 

which indeed generate Darboux transformations. The superhamiltonian and the two su- 
percharges form the so called SUSY algebra. Such approach appeared to be very useful 
in studying a variety of different kinds of quantum systems. A very fruitful approach 
is based on generalizations of one specific part of the SUSY algebra, namely the invari- 
ance of the Superhamiltonian in respect to supertransformations. For the SUSY partner 
Hamiltonians this invariance is expressed by intertwining relations. 

Below we list some specific classes of problems which were investigated in this way: 1) 
the class of shape-invariant potentials ^ with the algebraic construction of the spectrum; 
2) supercharges Q^ of higher orders (HSUSY) in derivatives 0-0 including irreducible 
2-nd order transformations and polynomial deformation of SUSY algebra; 3) inter- 
twining of Hamiltonians H12 with complex potentials by complex supercharges up to 



second order [|T0|; 4) construction of integrable 2-dim quantum and classical systems and 



corresponding symmetry operators |[Tl| ; 5) investigation of systems with matrix potentials 



In this paper we investigate some further generalizations of 1-dim SUSY intertwining 
relations starting from two hermitian Hamiltonians intertwined by complex supercharges 
of second order in derivatives (Section 2). A system of intertwining relations arises so that 
the same Hamiltonians are related at the same time by first and second order Darboux 
transformations. The compatibility of these constraints require the potentials to satisfy a 



stationary KdV equation |T3|. This formulation automatically allows (by elimination of 
one of the partner Hamiltonians) for a discussion of symmetry operators Ri of third order 
for each partner separately [Hi, Rj\ = 0. These symmetry operators which commute with 
the Hamiltonians were studied before in another approach |T^,[|I1]5 ^^^ the corresponding 
potentials were found. 

In Section 3 a further generalization of the previous relations for the symmetry op- 
erators is obtained by the ladder equation [H, a] = 2Xa. This approach reflects in part 
results obtained in |T6[ in their construction of dressing chains in connection with in- 



tegrable bi-hamiltonian systems and in [|1^ in the context of numerical illustrations of 



(cyclic) shape-invariance preserving standard shape invariance in each step. However 
our approach also leads to new findings concerning explicitly with the role of irreducible 
second order Darboux transformations 0. We will focus attention on the non-trivial de- 
pendence of the spectrum on the relations between particular parameters of the system, 
illustrating in detail the particular cases of second and third orders. We find for the third 
order different type of spectra involving up to three sequences of equispaced levels while 
the potentials are not strictly harmonic. Conditions for which one of these sequences is 
shrinked to one or two levels are investigated via the study of zero modes of the gener- 
alized creation operator. We stress that these zero modes are allowed only because the 
creation operators are here of higher order. 

In the Conclusion (Section 4) we emphasize the original aspects of our investigation 
in contradistinction with results in the literature. 

2. Complex supercharges, symmetry operators of third order 
and KdV equation 

In this Section we introduce supercharges with complex superpotentials intertwining 
two hermitian partner Hamiltonians (|l|), where: 

H, = -d' + Vi{xy, d^^. (2) 

ax 



Since complex supercharges of first order do not lead to a genuine extension we start from: 

Q+ = M++iq+; (3) 

M+ = d"^ -2f{x)d + b{xy, (4) 

g+ = -2f{x)d + b{x). (5) 

As a consequence of (JlD we obtain that f{x) = Const and one can rescale the first order 
g+ = — 2/g+ with g+ of standard structure: 

q^ = d + Wix); Wix)^-M. (6) 

2j 

In terms of the second order M"*" and the first order g"*" (|I|) amounts to a system of two 

intertwining relations for the same partner Hamiltonians: 

HiM+ = M+H2; (7) 

H,q+ = q+H2. (8) 

Imposing (J^) leads |^ to: 



,2,... , /"(^) /"(^) <* 



V,,W = T2m + /^W + ^-^-jj^ + a; (9) 

M.) ^ -/'(.) +m-^ + ^.j^ (10) 

with arbitrary real constants d and a. Let us remind that the sign of d in (^), (|iy) is 
crucial to characterize the so called reducible {d < 0) and irreducible {d > 0) second 
order transformation Q. 

Similarly, (H) gives the standard 

Fi,2(a;) =±W^'(x)+iy2^x), (11) 

where an additive constant is set equal to zero. 
Consistency requires that 

W{x) = -2f{x) +7; 7 = const (12) 



^We remind that reducible second order operator Af+ can be written as a product [d + Wi{x)){d 
W2{x)) where the real superpotentials Wi[x),W2{x) satisfy a gluing condition —W[{x) + Wl{x) ■ 
+W2{x) + W^ix) + y^ (see details in 0). 



and that f{x) satisfies a nonlinear second order differential equation which can be inte- 
grated to: 

f^(x)=4f\x)~8-ff{x)~4{a--f^)f{x)+Pf{x)~d; [5 = const. (13) 

Though ( pTSD can be solved for x = x{f), we will not make use of this solution. 

It is known [H], |]T^, that for arbitrary intertwining relation for hermitian Hamilto- 



nians one can construct symmetry operators of even order in derivatives which can be 
more or less trivial in 1-dim scalar case but not for 2-dim or matrix case. Since now 
we deal with two intertwining relations (0), (^), we can provide non-hermitian symmetry 
operators Ri = M"*"g~,i?2 = M~q^ of third order, to which of course one can add an 
arbitrary constant: 

HiM+q- = M+q-Hi; (14) 

H^M-q-^ = M-q+H2. (15) 

One finds that the hermitian part of Ri is trivially proportional to Hi but its antihermitian 
part Ai is more involved: 

A, ^ -d^ + (Qf-3f'-6^f-a + ^')-d+i-12f + 18^f + Qia-^')f + 6ff-3^f) (16) 

with similar expression for A2. 

A systematic investigation of symmetry operators of higher order in derivatives was 



performed for 1-dim Schrodinger operator in [14], following the basic strategy of writing an 



operator polynomial in derivatives and imposing that it commutes with the Hamiltonian. 
By this procedure one obtains that third order symmetry operators exist if the potential 



satisfies the stationary version of KdV equation |13 



V"{x) = 3V^{x) - ciV{x) + Co (17) 

where co,ci are real constants. One can check that symmetry operator (|16D with f{x) 
satisfying (|13[) can be rewritten as 

A^d' + i-lv^{x) + j)-d- Ivlix) (18) 
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which indeed coincides with the expression obtained in [Q for potentials Vi (x) for which 
(0) holds. 

In general, the solution of (|l^ can be expressed |jl3| in terms of periodic Jacobi elliptic 
function: 

V{x) = Kisn'^{K2x;k) + K3, (19) 



where the constants Kj are related []T^ to the constants q of (|T7p. We focus attention now 
to the interplay between existence of a symmetry operator and the "non-degeneracy" 
of normalizable solutions in 1-dim Quantum Mechanics. For this purpose we restrict 
ourselves to the period k = 1, which corresponds to sn{K,2x;k = 1) = th{K,2x). It is 
straightforward to check explicitly that the only bound state of the associated potential 
Vi{x) is an eigenstate of Ai. More generally, as suggested by a fundamental theorem in 
Mechanics, one can convince oneself that the operator Af is a third order polynomial with 
constant coefficients in Hi for Vi{x) satisfying ([T7|) . 

As a final comment let us mention how one can extend the third order results to 
higher orders, for example to the fifth order symmetry operators. Imposing directly the 
commutation between the fifth order symmetry operator and the Hamiltonian, one obtains 
as a consistency relation that V{x) has to satisfy a higher KdV equation |]13| and that 
this symmetry operator can be written as a product HAi with Ai given in (|TB|). One can 
notice that if V{x) satisfies ([T7|) it will also automatically satisfy the higher KdV (but not 
vice versa). 

3. Intertwining with shift and higher order shape-invariance 

For illustration of the construction let us start by modifying equations (|^, (H): 

HiM+ = M+H2; (20) 

Hiq+ = q+{H2 + 2X), (21) 

introducing thereby a shift by the positive constant 2A. Eliminating one partner Hamil- 
tonian (for example, H2) and denoting product operators a~^ = q^M~ and a~ = M^q" 
we obtain: 

Hia+ = a+{Hi + 2\). (22) 



We will call such Hamiltonian as "third order shape- invariant" . One can also work the 
other way around and start from (|22|) with an operator a+ of third order represented as 
a product q^M~ and Hi represented by Hi = q^q^ as is familiar in SUSY Quantum 
Mechanics. After one can introduce an auxiliary Hamiltonian iJ2 + 2A = q^q^ to obtain 

Eq. (|2^) is a ladder equation where a"*" plays the role of generalized creation operator 
which provides an excitation energy of 2A. In order to study the spectrum it is crucial 
to study zero modes of a~ and a"*". The former describe the lowest lying levels of the 
system, and one has to apply recursively the operator a^ to them in order to generate 
the excitation spectrum. The energies of the zero modes can be obtained by imposing the 
vanishing of their norm, which involves the average of the operator product a^a~ . This 
product can be easily evaluated algebraically because 

a+a- = q+M-M+q- = q+{{H2 - a)^ + d)q- = Hi{{Hi -a- 2Xy + d). (23) 

Contrary to the standard harmonic oscillator, one has the possibility to have also zero 
modes of the operator a^ which correspond to a possible truncation of the tower of excited 
levels. Arguing as before, the relevant operator product is in this case 

a-a^ = {Hi + 2A) {{Hi - a)^ + d) (24) 

to be averaged over the excited states. 

This construction can be generalized to higher orders taking into account that on gen- 
eral grounds n— th order operators a^ can always be constructed in terms of products 
of q and M. In the general higher order case, (p2D gives a connection between Hi and Hi 
plus shift, which is the simplest realization of the notion of shape-invariance |^. Though 
some properties of the spectrum, like zero modes of a^, will depend on the explicit prod- 
uct structure of a^, the excitation spectrum can be mainly obtained algebraically. So 
it is possible to study the consequences of ( ^2] ) without taking into account the specific 
definition of the operators a^. 

In order to provide additional arguments to the interpretation of Eq. (p2|) as a gener- 
alization of shape invariance, let us start to notice that if a^ would be of first order. Hi 
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would be identified with harmonic oscillator. If a"^ is of second order, we reobtain the 
singular harmonic oscillator potential which is also shape-invariant (see next Subsection). 
We stress that if a^ would be an operator of second order, it is known that it is not 
always possible to write it as a product of two first operators with real superpotentials: 
such case is referred as irreducible. 

3.1. Intertwining with shift and second order shape- invar iance 

In this Subsection to keep notations the same we consider (PO]), ( PT| ) with M"*" tem- 
porarily of first order. We can thus explore in this Subsection if even a first order in- 
tertwining operator M"*" can now lead to nontrivial consequences. It is easy to convince 
oneself that this is indeed so. The superpotentials which solve (PDj), ( PT] ) are superposi- 



tions of a term growing as x and a singular term like 1/x leading to two singular potentials 
g Vi{x) and Viix) as foUowsQ: 



p{p - 1) X'x' p{p + l) X'x' 

~^ ^ ^~5 V2(x) = \ - 

x^ 4 x2 4 



^.M^i-^ + n^; ^,M . i-^ + 11^ _ A. (25) 



These potentials are shape-invariant in the standard sense f^ and belong to the class of 
algebraically solvable models, because V2{x; p. A) = Vi{x] p + 1, A) — A. 

It is instructive to display the algebraic properties of these systems. In terms of the 
product operators a'^ = q^M~ and a~ = M^q^ (we repeat that in this Subsection a^ 
are of second order) one obtains an algebra suggestive of a generalization of the standard 
harmonic oscillator algebra: 

[i/i, a+] = 2Aa+; [a+, a'] = -AXHi + const. (26) 

Similar results hold for H2. 

It is now tempting to consider, as explained before, Eqs.(PB|) by themselves including 

the case where the operators a^ are irreducible. We know from 0] with H2 = Hi + 2A 
^On the half hne for p = I + 1 and suitable choice of boundary conditions for wave functions 
potential Vi can be interpreted as a radial harmonic oscillator in partial wave I. 



and setting^ a = in (Q) that 

a-^a- = Hi + d] a'a^ = {Hi + 2Xy + d. (27) 

Solving (|2^) we obtain from (||) that /(x) = |Ax and 

Vl(^) = ^-(4 + ^)^-A. (28) 



Evaluating the matrix element of the first of (|27[) , we get the equation for the 0-modes 
of a^: El 2 = ±a/— rf. For a positive value of d (irreducible case) there are no solutions 
and physically this is clear because the potential ( pSj ) describes a collapse for d > 0. 
For d < the singular part of the potential is less pathological and one may have two 
zero modes for suitable choices of parameters. The eigenfunctions of the zero modes of 
a~ which are also eigenfunctions of Hi can be found explicitly by replacing d^ in a" by 
(Vi(x) —E^^^). The asymptotic behavior of the corresponding eigenfunctions is compatible 
with the normalizability. 

We expect that for self-adjoint extensions of Hamiltonian Hi both operators a~^a~ and 
a~a^ are nonnegative on physical states. There exist two zero modes of a~ if \/—d < X. 
Acting with a+ on these two zero modes, one can construct two sequences of levels with 
internal spacing given by 2A. In the limiting case d = the two sequences coincide. 
When ^/—d > A, one obtains only one zero mode with the energy E^^^ = \/ —d and, 
therefore, only one sequence of equispaced levels. The existence of the second zero mode 
E^^' = —\J—d in this case is incompatible with positivity of the second operator in 
Eq. (P7[) . For \J —d = A this second zero mode remains non-normalizable. 

3.2. Intertwining with shift and third order shape- invar iance. 

From now on we will consider the intertwining operator M^ to be a reducible or 

irreducible operator of second order (^ but q^ still of first order @. The solutions of 

(pop are unchanged in respect to Eqs.(|), (p!OD of Section 2 but Eq.(|2T|) implies a shift 

for the potential V2{x) in ([TT]). Thus the consistency equations for W{x) and f{x) are 
■^This choice of the energy scale means M that the eigenvalues of Hi are bounded from below by 
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modified. Eq.([T2|) becomes: 

W{x) = W3{x) = -2f{x) - Xx, (29) 

where an additional integration constant can be ignored because of a sliift of x, whicli 
fixes tlie origin of coordinate^. The potential Vi{x) can be written from (JTlJ 



as: 



Vi{x) = -2f'{x) + Af\x) + AXxfix) + AV - A, (30) 

with /(x) satisfying: 

r = Qrl + 6/'(^) + 8\xfix) + 2(AV - (A + a)) fix) + -f--. (31) 

The equation ( ^Tj) can be transformed by the substitution f{x) = l/2V\g{y); y = VXx 
to the Painleve-IV equation [|T9[: 



where 



n''^(v] 3 h 



(1 = 1 + ^; ''=-^- (33) 



This equation has been studied intensively in the last years [EDI and in the following we 



will be mainly interested in asymptotic properties of its solutions which will determine 
the asymptotics of potentials (|30D and the normalizability of eigenfunctions. 

Concerning the algebra, the only modification in respect to (^) in the previous Sub- 
section is given by: 

[a+, a-] = -2\{3H^ - (4a + 2\)Hi + a"^ + d) (34) 

and similarly for H2. When the index will not appear explicitly we intend to refer to the 
index 1. 

We show how one can derive the spectrum from (^) and (|3^ if normalizable zero 
modes of the annihilation operator a~ exist. We stress that this algebraic method is very 



"'in Eq.(p^) there was no such fixing, so the hmit A ^ cannot be taken to obtain the results of 
Section 2. 
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powerful now since the explicit form of the potential is known only in terms of Painleve 
transcendents. 

The equation for zero modes of a~ reads: 

a-^i°^ = M+g-^^°^ = (35) 

where k labels the normalizable solutions. 

As explained at the beginning of Section 3, Eq.(^), the algebraic equation for eigen- 
values El is: 

^(0) . [(^(0) _ ^ _ 2A)2 + d]=0, (36) 

and has at most three real solutions. 

For the irreducible case (d > 0) only one zero mode Eq = exists^. But it is necessary 
to check indeed the normalizability of corresponding zero mode which can be written (|35D, 
(PPP, (^) in terms of superpotentials: 

^f,°^(x) = exp{rW{z)dz) = exp{--Xx^ - 2 r f{z)dz) (37) 

where f{z) satisfies (0) and must have a suitable asymptotic behavior in order to make 
(P7|) normalizable. It has been found [^ that for 6 < there are two possible asymp- 
totic behavior of g{y) at ±oo : —2y/3 and —2y. It has been demonstrated [^ that 



one can realize the matching between the same asymptotic value at +oo and — oo. The 
asymptotic values which are compatible with normalizability of ( PTf ) is g{y) ~ —2y/3, 
or equivalently, f{x) ~ —Ax/3. We notice that the subleading term in f{x) have oscilla- 
tions ~ cos(Aa;^/v3 + o(x)) which show up in the potential (|30|) as a modulation of the 
harmonic term. 

From [H, a"*"] = 2\a~^ one can deduce that E^"'^ = 2\n. The related potential (|30|) for 
the increasing asymptotic behavior of f{x) has asymptotics X'^x'^/9, which differs from 
conventional A^x^. We remark that the latter system does not belong to the families of 
potentials isospectral to harmonic oscillator pl,[^ which have the same asymptotics but 
differ from harmonic oscillator for finite x. 



'^Tn -fViTC riooo it7£1 -n/^i-fi/^o fVio+ cinr'Ci A/T aT/ * ' 



In this case we notice that since M ^g 7^ also H2 has the ground state with zero eigenvalue. 
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Evaluating the norm of the excited state '^^"'~^^\x) = a"'"\l/(")(x), one finds after alge- 
braic manipulations explained in the beginning of Section 3 that for the irreducible case it 
cannot vanish. So, the spectrum is not bounded from above. Analytic calculation support 
clearly this result because of the asymptotic oscillator-like growth. 

For the reducible case M"*" with d < is: 

M+ = {d + Wi{x)){d + W2{x)) (38) 

with real superpotentials 0] 

W^M = -f{x)±l^jj^. (39) 

Eq.(p6D has the solution Eq = as before and two additional solutions: 

e'l^ = a + 2A - V^ (40) 

^f = a + 2A + V^ (41) 

which correspond to eigenvalues oi H = q^q^ provided (0), (|41|) are non-negative and 
that associated eigenfunctions \E'^ (x) are normalizable. Acting on the three zero modes 
W^ {x) by the creation operator a"*", one creates the excited states of the system or- 
ganized in three sequences. According to Eq . (|2^) , within each sequence the levels are 
2 A— equidistant. The non-negativity of the norm ||a^\l'(x)|p for all physical states '^{x) 
leads to additional necessary conditions for the spectrum^. 

The eigenfunctions \E'^ (x), annihilated by a~ = M^q" = {d + Wi{x)) ■ {d + W2{x)) ■ 
{—d + Ws{x)) with Ws{x) defined in (p9|), can be calculated exphcitly: 



^{,°^(x) = exp{ W3{x')dx'y, 

^5')(x) = {W2{x)-W3{x))exp{- rW2{x)dxy, (42) 

*L°^(a;) = [2V^+{W2{x)-W3{x)){Wi{x) + W2{x))]exp{- rWi{x')dx'). 

Their normalizability, similarly to \l'o (x) for the irreducible case, depends on the choice 

for the asymptotics of f{x) and particular values of parameters a, A and d. In reducible 
^From the algebra (p3) - (EJ) one can prove that these conditions are sufficient to ensure the non- 
negativity of the norm of any state created by a polynomial of a^ : ||P(a+, a^)^(a;)|p > 0. 
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case b > four possible asymptotics of solutions of Painleve-IV equation can appear: 
giy) ~ -2y/3, -2y, and ±Vb/2y. 

On the other hand, the operator a^ can also have zero modes which can be found 
explicitly: 



\l'i(x) = exp{ Wi{x)dx); 

^2(x) = {Wi{x) + W2{x))exp{rW2{x')dx'); (43) 

^3(2;) = (e^^^ + {Wi{x) + W2{x)){W2{x)-W3{x))yxp{- I^Ws{x')dx'). 

These eigenfunctions are solutions of the Schrodinger equation with eigenvalues: 

Ei = a- V^; E2 = a + V^; E^ = -2A. (44) 

The consistency of the energy spectrum, i.e. the non-negativity of operators a'^a~ and 
a~a~^, rules out the possibility to have a normalizable zero mode with negative energy E^ 
for nonsingular superpotentials. Concerning the other zero modes, one can conclude only 
that the total number of zero modes of both operators a~ and a~^ cannot exceed three, 
which follows straightforwardly from the conflicting asymptotics (H^) and 



As the spectrum of Hamiltonian is bounded from below and the operator a"*" raises its 
levels, zero modes of a"*" represent an obstruction to build an infinite sequence of levels. 
Only one zero mode of a"*" may exist together with two zero modes of a~. In this case 
the spectrum consists of one infinite sequence of levels and a finite band of levels, both 
generated by operators (a^)". 

Let us proceed to a more detailed analysis of different spectrum patterns. 

a) Three normalizable zero modes (^) of a~ and, respectively, three equidistant 
sequences of levels may arise only for A > \/—d if \^—d < a < 2X — \/—d or yj —d—2X < 
a < —\/—d. The corresponding solutions of Painleve-IV equation must have the leading 
asymptotics f{x) ~ —Ax/3 with subleading oscillations cos(Ax^/v^ + o{x)). 

b) Two normalizable zero modes of a~ (and two sequences of levels) may exist for 

different solutions of Painleve-IV equation Q. Namely, for A > \/~d the solution with 

^We discuss here only the case of equal asymptotics of f{x) at ±cxd. However we are aware of existence 

of solutions with the different asymptotics, which may lead to an additional variety of potentials with one 
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asymptotics f{x) ~ —Ax/3 provides the fall off of '^_'{x), \E'Y (^) if a > 2A 



(0), 



,{0)/ 



-d 



M 



(0), 



and of ^r(a;), ^r(a;) if -2A 



-d < a < — 2A + \J —d. For A > \J —d the solution 



with asymptotics f{x) ~ —\/—d/2\x generates two sequences of levels starting from 
^{,°^(x), ¥^\x) if 2An + v^^ < a < 2A(n + 1) - V^] n = 0,1,2... For A < y/^ 
two sequences (from \1'q (x) and \1/+ (x)) are generated for two possible asymptotics 
/(x) ~ v^— ^/2Ax and /(x) ~ —Ax/3 if — 2A — ^/—d < a < —\J—d. 

c) For a < — 2A — ^/—d and arbitrary positive A one sequence of equidistant levels will 
be realized^ with ground state \l/o (x) and one of three asymptotics: /(x) ~ —Ax/3, ~ 
±^/—d/2Xx. For /(x) ~ —Ax/3 only one sequence of levels will be realized also for 
A = \/—d, a = ±y/—d (starting from \1/q (x)), for A = \/~d, a > —\/~d (starting from 
\1'_ (x)) and for A < \/—d, a > — 2A — \/—d (starting from \1'^ (x)). For asymptotic 
behavior /(x) ~ —^/—d/2Xx and A > ^/—d the spectrum consists of one sequence of 
levels if \a — 2Xn\ < V—d, n = 0,1,2,... (starts from \l/_ (x)) or if a = \/—d (starts 
from \E'o )• Last, for /(x) ~ \/—d/2\x one obtain one sequence with ground state \l/\_ if 
A < \^—d and a > —\/—d or if A > \/ —d and \(x — 2\n\ < \/—d, n = 0,1, 2, ... 

For specific values of parameter a the existence of normalizable zero mode of a"*" can 
truncate one of two sequences ( with ground state \l/o (x) ) in item b). 

d) One additional singlet state PI satisfies the equation 



a+^ 



(0), 



X 



a"\l/ 



(0), 



X 



0. 



For A > v^^ (A < y—d) it occurs when a = =i= 
entails the equation: 

/'(x) = -2f\x) - 2Ax/(x) T v^ 



d and ^{,°^(x) 



^2,1(2;) which 



(45) 



with asymptotic behavior /(x) ~ ^\/—d/2\x. The spacing between the two ground states 
L°^(x) is: AE = 2AT2y^. 



^^q\x) and^ 



or two sequences of energy levels. The admissible range of parameters and zero mode space is derived 
then by intersection of those ones for the particular asymptotics described below. 

^Let us remark that particular value a = — 2A — %/— d, i.e. E\^ = 0,E_ < 0, just corresponds to 
the case when Painleve-IV equation (p2) has a class of particular solutions which coincide with solutions 



of the Riccati equation g'{y) — g^{y) + 2yg{y) + \fh. Substituting this Riccati equation into Eq.(30) one 
finds that potential is the pure harmonic oscillator: V\{x) = }?x^ — A. 
^ Its existence has been mentioned in [ESl 
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e) The doublet representation {Wq (x), a^^o (x)) of the spectrum generating alge- 
bra (P^), ( P^ ) is built on solutions of the equation: 

It may hold when a = 2\ + \J ~d for arbitrary positive value of A and when a = 2A — \/ ~d 
for A > \/—d. It is equivalent to cl^'^q (x) = \l/i 2(a;), which is satisfied when f(x) obeys 
the following equation: 

8Xfix)if{x) + 2 fix) + 2Xxfix) ~2Xt V^) = (fix) + 2f{x) + 



2\xf{x) T v^) 



{fix) + 2f{x) + 2\xf{x) -2\T V^) 



(-fix) + 2f{x) + 2\xf{x) - 2A T v^^) - 4A(A ± v^) 



(46) 



One can show that all solutions of this equation fulfill the Painleve-IV equation (|3lD. These 
solutions have the asymptotics f{x) ~ ±-\/— d/2Ax and cannot have any (pole) singularity. 
The spectrum consists of a doublet (0, 2A) and infinite sequence E^ = ±2\/—d + 2{n + 
2)A, n = 0,1,2... 

It may seem to be possible to obtain a higher representation of spectrum generating 
algebra with three or more levels in the band. However for higher order Darboux trans- 
formations it is known that the zero mode of the third order operator a^ cannot have 
more than two nodes, i.e. it cannot be higher in energy than the second excited level 
of system. From our analysis it follows that states forming the finite representation lie 
below the first state of the infinite sequence of levels, which is obviously a zero mode of 
a~ . Therefore this finite representation may include only the ground state and, possibly, 
the first excitation. For shape-invariance of A^-th order the maximal dimension of finite 
representation may be {N — 1). 

4. Conclusions. 

As conclusive remarks let us now emphasize the novel features of our approach. We 
would like to point out what in our opinion is new in respect to |l^. First, we introduced 
in the dressing chain (in their language) irreducible transformations which correspond 
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to quantum mechanical systems with only one sequence of levels. Second, we discov- 
ered that in suitable circumstances one sequence can be truncated because the higher 
order creation operators a^ can have zero modes. Third, we have given attention to the 
asymptotic properties of the Painleve-IV equation (^) in connection with discussion of 
the normalizability of eigenf unctions. In particular, for the construction of potentials 
which do not have subleading oscillatory behavior, we revealed the importance of the 
special asymptotics g{y) ~ Vb/2y which can be implemented [^ at both ±00 or only 
alternatively at +00 or at —00. We have found that a particular class of solutions of a 



Painleve-IV equation can be obtained as solutions of another, first order, equation (^Bf) 
which does not admit singularities. 



Concerning [|I7|, we differ because implementation of higher order shape-invariance 
in our scheme does not require to preserve standard shape-invariance of first order step 
by step. The spectra that we can generate appear to be more rich: possibility of having 
three, two, one sequences; possible truncation of one sequence; non-oscillatory asymptotic 
behavior of the potential. 

In respect to our approach gives more importance to the distinction between re- 
ducible and irreducible second order transformation and to potentials which generate 
sequences of levels spaced by 2A but do not in general have the asymptotics behavior 
AV. 
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